EUROPHYSICS LETTERS 
Europhys. Lett., (), pp. () 



Wannier threshold law for two electron escape in the pres- 
ence of an external electric field 

Bruno Eckhardt 1 and Krzysztof Sacra 1,2 

1 Fachbereich Physik, Philipps-Universitdt Marburg, D-35032 Marburg, Germany 

2 Instytut Fizyki im. Mariana Smoluchowskiego, Uniwersytet Jagiellonski, 
ul. Reymonta 4, PL-30-059 Krakow, Poland 

(received ; accepted ) 

PACS. 34. 80. Dp- Atomic excitation and ionization by electron impact. 
PACS. 32.60. +i - Zeeman and Stark effects . 

PACS. 31.25.Jf - Electron correlation calculations for atoms and ions: excited states. 



Abstract. - We consider double ionization of atoms or ions by electron impact in the presence 
of a static electric field. As in Wanniers analysis of the analogous situation without external 
field the dynamics near threshold is dominated by a saddle. With a field the saddle lies in a 
subspace of symmetrically escaping electrons. Near threshold the classical cross section scales 
with excess energy E like a ~ E a , where the exponent a can be determined from the stability 
of the saddle and does not depend on the field strength. For example, if the remaining ion has 
charge Z = 2, the exponent is 1.292, significantly different from the 1.056 without the field. 



The threshold behavior of double electron escape from an atom or ion was first tackled 
by Wannier in a famous paper in 1953 [El 0, |1. On the basis of the classical dynamics of 
two electrons in an atom he concluded that on account of the electron repulsion the two 
escaping electrons are correlated and that the cross section increases with excess energy E like 
o~(E) oc E a , with a non- integer exponent a that depends on the charge of the remaining ion. 
Since the original paper many fine details of the process have been elucidated in both theory 
and experiment (for recent reviews see [5], ||). 

Electron-electron correlations are also important in the time-dependent process of double 
ionization in strong laser pulses. Measurements of the ion and electron momenta show that 
in non-sequential double ionization the escaping electrons prefer symmetry related motions 
[Q> h> hI' The processes that are important for this double ionization are a matter of debate, 
but when discussed within the rescattering model 10, nj similarities to the Wannier problem 
show up Jl2], [ll| Q . In the rescattering model one electron is temporarily ionized by tunneling, 
but driven back to the atom when the field changes the phase. During this rescattering event a 
highly excited two electron complex close to the nucleus is formed which then decays towards 
double ionization. For the field strengths where the characteristics of correlated electron 
escape are observed it seems crucial that the external field does not vanish when the decay 
takes place; otherwise the electrons do not have enough energy for double ionization [|l2| p"3| . 
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In our previous presentations of the process in a time-dependent field we argued that the 
electron motion is fast compared to the field oscillations so that an adiabatic analysis can be 
applied. But it is possible to discuss the process also in the presence of a static field, where 
no adiabatic assumption is needed, and to push the analogy to the Wannier problem further 
by deriving the threshold laws for non-sequential double ionization in the presence of a static 
field. This is our aim here. 

Wanniers analysis divides into two parts: the identification of the configuration that leads 
to double ionization at the threshold and the determination of the threshold law from the 
stability exponents of the fixed point. In Wanniers case the threshold configuration had two 
electrons escaping symmetrically on opposite sides of the nucleus, thus minimizing electron 
repulsion. The presence of the field introduces a preference for motion in the direction of the 
field gradient. At the threshold for the process the energy is equally distributed between the 
two electrons. Furthermore, the distance of the electrons to the nucleus should be the same 
for otherwise any difference would be amplified by electron repulsion and the electrons would 
not escape simultaneously. Therefore, in the presence of an electric field the configuration 
that corresponds to Wanniers is one where the electrons escape along trajectories which are 
reflection symmetric with respect to the field axis |12, O] . 




Fig. 1. - Potential energy in the symmetric subspace, see Eq. (|lj), for Z — 2 and F = 1. When one 
decreases the amplitude of the field the saddle moves along the dashed line. 



Due to the rotational symmetry around the field axis the component of angular momentum 
along the field axis is conserved and for the configuration with lowest threshold this component 
vanishes. We can therefore analyze the process in a subspace of zero angular momentum with 
symmetrically placed electrons. In this subspace we introduce cylindrical coordinates with 
z\ = Z2 = z the distance along the field direction and p\ = pi = r the transverse distance. 
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The electrons are on opposite sides of the core with respect to the field axis, so that ipi = (f> 
and if2 = <t> + 7r. The momenta are p Zl — p Z2 = p z /2, p Pl = p P2 = p r /2 (the factor 2 in the 
momenta is related to the proper canonical transformation and restriction of the dynamics to 
the symmetry subspace ]l3[]) and p Vl = p V2 — 0. With these coordinates for the symmetric 
subspace the Hamiltonian becomes, in atomic units, 

pi + pi 2Z 1 

H = Fr / z - , + — -2Fz, 1 

4 Vr 2 + z 2 2r w 

where P 1 is electric field strength and Z the charge of the ion after removal of the two electrons 
Equipotential curves for Z = 2 are shown in Fig. [|. Clearly noticeable is the saddle 

at 

- _(2«"1) 1/4 Ul _ (2a -1)3/4 



of energy 7 a = -2(2a - 1) 3/4 ylH where 

a=(2Z 2 ) 1 / 3 . (3) 

The ratio r a /|z a | depends on the charge of the nucleus but not on the field strength. Thus, 
when the field strength varies the saddle moves along a line that forms a fixed angle with the 
field axis, as indicated in Fig. [l]. 

The potential shows that in the symmetric subspace double ionization requires crossing of 
the saddle. In the full phase space this saddle acquires an additional unstable eigenvalue, so 
that the symmetric subspace is an unstable subset of the full phase space. Nevertheless, it 
controls double ionization events since trajectories leading to non-sequential double ionization 
must pass close to the saddle in the symmetric subspace: the instability of the subspace is 
connected with an increasing asymmetry in position, momentum and energy between the two 
electrons. Near the threshold for double ionization this results in trajectories being pushed 
towards single rather than double ionization. 

The stability analysis of the saddle in full configuration space gives one neutral direction, 
connected with an overall separable rotation around the field axis, three stable directions and 
two unstable ones: one unstable direction is the 'reaction coordinate' |ll| |l6) across the saddle, 
clearly visible in the potential in Fig. [j], and the other corresponds to motion away from this 
subspace. In the vicinity of the saddle the potential can be expanded to second order in the 
deviations from the saddle and the Hamiltonian can thus be approximated as 

E „v. + A-& + A-!& + ±(!L + i &y (4) 

»=i v 7 

where x denotes the reaction coordinate in the symmetric subspace, y the unstable mode away 
from the subspace and Ui are the stable modes with frequencies u>i. For the case of zero total 
angular momentum projection on the field axis we consider here there is no contribution to 
([|) from the neutral mode. For later reference we note the eigenvalues of the two unstable 
directions: for the reaction coordinate it is 



" ~ (2a -1)V4 < 



and for the motion away from the symmetric subspace it is 



2 _ y/32a - 28 + 6/a + 2x72a - 1 f3/2 
V ~ (2a -1)5/4 * ■ W 
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The escape from the subspace y is such that p\ = r s + c r y and z± — z s + c z y increase but 
pi = r s — c r y and z-i = z s — c z y decrease. The constants c r and c z have the same sign and 
determine the direction of the unstable mode in the configuration space (see jl3| for further 
details). 

For the energy equal to the saddle energy only a trajectory living in the symmetric subspace 
can lead to simultaneous double escape. This reduces the dimensionality of the process below 
that of the energy shell and consequently the cross section vanishes. For higher energy some de- 
viations from the symmetric configuration are allowed, bringing the subset to non-zero measure 
on the energy shell. The energy dependence of the cross section near threshold is determined 
by a competition between the two unstable modes fl7|| . This is in close correspondence to 
the motion on the Wannier ridge where only trajectories sufficiently close to the symmetric 
configuration can leave the Coulomb zone and ionize [[I). 

All trajectories leading to the simultaneous escape of the electrons have to pass near 
the saddle. Thus, to estimate the cross section close to the threshold we may employ the 
Hamiltonian of the system in the harmonic approximation, Eq. (^). The cross section can be 
calculated from the phase space flux associated with the trajectories that lead to double 

ionization. With V the projector onto these trajectories, the microcanonical phase space flux 
j a at energy E that crosses the saddle can be calculated at the surface x — with the phase 
space velocity projected onto the normal of that surface, 

3 

jo = Px pS(x )S(E - H)V dp Xa dx dp VQ dy dp Ui0 du i0 

i—l 

= / pVdpy o dy Y[dp Ul0 du iQ , (7) 

i—l 

where p = p(p X ai x OiPyoiyo,Pwi T w io) is the initial phase space density. 

We are interested in the dependence of the cross section on the energy e = E — V s above the 
saddle. The integration limits of the stable degrees of freedom are arbitrary but finite, so that 
they cannot go to infinity, as required for the escape from the reaction zone towards ionization 
[Q. The stable directions thus do not contribute to the energy dependence, just as in Wanniers 
example [0, [l8| ^7) . Thus the only degree of freedom with a critical energy dependence in (Q) 
is the unstable mode associated with yo and p Va . Assuming the initial phase space distribution 
to be approximately uniform and energy independent Q we get jo-(E) cx J Vdp yo dyo. 

Taking the initial reaction coordinate on one side of the saddle, e.g. xq < with p Xo > 0, 
the initial conditions of the unstable mode must be chosen so that after crossing the saddle 
both electrons escape. Since x and y are center of mass and relative coordinates, respectively, 
double escape requires that both x + y and x — y simultaneously go to infinity. Thus, the 
projection operator V selects those orbits that at some distance x after crossing the saddle 
satisfy 

\y\<\x\. (8) 

From this it follows that initially only a small amount of energy e — E — V s can be put in 
the y mode because the Lyapunov exponent for the desymmetrization is larger than that for 
motion along the reaction coordinate, v > p. Hence for small e, the initial momentum p Xo can 
be approximated as 



Px = J^e + p 2 xl ~p 2 ya +v 2 yl » p\x \ + - — r. (9) 
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Expressing y in terms of the initial conditions gives, for large time, 



.,!«)+—)( r i ) v/ " ,;:Mii 

2 V v J \x +p xo /n; 



1 / ( 2x 



Substituting (|l0|) and @ into (||) then results in 



?/o + — 
v 



< const • e^, (11) 



which is precisely the restriction we need. Changing to canonical variables y' = y + p y /v 
and p y i = Py/2 — uy/2, one finds that p y >{t) — p y iQexp(—ut) while y'(t) = y' exp(i>t). The 
inequality (Fyl) is thus the restriction on the unstable direction in the phase space. The initial 
condition along the stable direction, p y 'o, can be arbitrary - it is, of course, finite due to 
the requirement that ionizing trajectories must emerge from the reaction zone [|lj. The cross 
section a(E) is proportional to the flux j a , so that we find for the threshold behaviour the law 



a(E) cx j„ cx J V dy' a cx (E - V s ) a , (12) 



with the exponent a given by 



2 _ v 2 _ v/32a - 28 + 6/a + 2y/2a - 1 
a ~ ~ ^/SOa- 49 + 12/a - ^J2a - 1 ' 

where a is defined in Eq. (|). In Fig. | we show a as a function of the nuclear charge together 
with the Wannier exponent; the values are also listed in Table I. For increasing charge Z the 
exponent a decreases and approaches the limiting value y/3/2 ~ 1.225. 

The exponent a does not depend on the field amplitude but only on the charge of the nucleus. 
This raises the question how the classical Wannier theory can be recovered for vanishing field 
strength F. The answer can be given easily by appealing to Wanniers analysis: he introduced 
a Coulomb dominated zone inside which potential energy exceeds kinetic energy. The outer 
boundary of this zone is then energy dependent and increases without bound as the total 
energy approaches zero. As long as the saddle induced by the field is inside this radius we 
expect the exponent given above (|l3|), but if it is outside we expect the classical Wannier 
exponent. In the limit of vanishing field strength but for fixed excess energy the saddle moves 
outside the Coulomb zone and the classical exponents are restored. This transition might be 
accessible experimentally. 

The classical Wannier theory for double ionization without external field gives the same 
threshold law as do semiclassical and quantum calculations |B This may be traced 

to the fact that a remarkable scaling law relates the limits of energy approaching threshold 
and of Planck's constant becoming small |Q, |J. In the present case the saddle is at a finite 
distance and can be overcome by tunneling. This will modify the threshold law very close to 
the classical saddle, in an energy interval of about hp,, with p, the frequency for motion along 
the reaction coordinate (eq. (5)). In the semiclassical limit and for weak fields this can be 
made arbitrarily small, so that the algebraic behaviour should become accessible. 

The most direct experimental study of the proposed threshold behavior would be doublc- 
photoionization || in the presence of an external field. The photon resolution of 0.1 eV 
achieved in recent experiments (|l0[ | should enable detection of deviations from Wannier theory 
e.g. in double-photoionization of He atoms in a static field of 30 kV/cm, where the saddle 
energy is —0.3 eV. 

Pulsed lasers provide stronger fields but add a time-dependence. However, the laser field 
can be considered as stationary if the time needed to cross the barrier is a small fraction of 



6 



EUROPHYSICS LETTERS 



1.4 




2 4 6 8 10 

Z 



Fig. 2. - Exponents in the threshold law for different charges of the nucleus. Full circles correspond to 
Eq. ( J13| ) in the presence of an external field while open symbols give the Wannier prediction without 
field |||. 



the field cycle, only. Double ionization could be triggered by a crossed electron beam. The 
energy of the saddle and of the electronic beam change with the field phase but both are well 
defined. The field intensity can then be adjusted so that the excess energy reaches the energy 
of the saddle only at specific moments in time. Such an experiment would also provide a very 
interesting direct test for the rescattering model for double ionization in strong fields Jl0|, [llj . 

Finally we should mention that the correlated escape used to derive the threshold law is not 
the only pathway to double ionization. Trajectory studies ]l2| , |l3| show that initial conditions 
started near the saddle that do not lead to immediate double ionization can lead to sequential 
ionization in that one electron escapes immediately but the other escapes after a return to 
the nucleus. However, in such a sequential escape the momenta of the second electron are 
not correlated to the ones of the first electron and this can perhaps be used to distinguish 



correlated from sequential escape i3 



To conclude we have presented the analysis of the threshold law for the double ionization 
in electron impact in the presence of a static electric field. Simultaneous ionization of the 
electrons then proceeds via a correlated crossing of the saddle induced by the external field. 
The cross section behavior close to the threshold is algebraic in excess energy with an exponent 
determined from the ratio of the positive Lyapunov exponents of the unstable modes of the 
saddle. We hope that the generalization of the Wannier theory to the case with an additional 
external field presented here will stimulate experimental and further theoretical investigations. 

Financial support by the Alexander von Humboldt Foundation is gratefully acknowledged. 
The work is also a part of KBN project No. 5 P03B 088 21 (K.S.). 
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